AD”A252  709 


4 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


COMPUTER  SIMULATION  OF  ACOUSTIC  FLUCTUATIONS 
DUE  TO  FINESCALE  TEMPERATURE  PERTURBATIONS 
MEASURED  BY  THERMISTOR  CHAIN 

by 

Gregory  D.  Crabtree 
March  1992 

Thesis  Advisor:  Pecheng  Chu 

Approved  for  public  release:  distribution  is  unlimited, 


92-18176 

■niiiiii 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


:i  s  sci  4*1 1 1 : 1 1 u  , 


REPORT  DOCUMENTATION  PAGE 


UNCLASSIFIED 


Approved  for  public  release; 
distribution  is  unlimited 


6c.  ADDRESS  (City,  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


Naval  Postgraduate  School 


7b.  ADDRESS  tCt^,  State,  and  ZIP  Code) 

Monterey,  CA  93943-5000 


f ^  fj  il  I  ^  J  ^  I  ill  ^  M  ^  : 1 1  ^  re 


I  f  ^  .T I  I 


ORGANIZATION 

(if  applicable) 

8c.  ADDRESS  (City,  State,  and  ZIP  Code) 

ELEMENT  NO.  NO. 


ACCESSION  NO. 


1 1  TITLE  (Include  Security  Classification) 

COMPUTER  SIMULATION  OF  ACOUSTIC  FLUCTUATIONS  DUE  TO  FINESCALE  TEMPERATURE  PERTURBATIONS  MEASURED  BY 
THERMISTOR  CHAIN  (U) 


1 

n*=iijflaMy,i:t;irj;va;MrAiw.iir 


14.  DATE  OF  REPORT  (Year,  Month,  Day) 

March  1992 


FROM  05/89  TO  03/90 


e  views  expressed  m  this  thesis  are  those  ot  the  author  and  do  not  reflect  the  ottici 
policy  or  position  of  the  Department  of  Defense  or  the  United  States  Government 


C06AT)  CODES 


GROUP 


SUB-GROUP 


18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

Acoustic  fluctuations,  finescale,  surface  duct  temperature  perturbations, 
parabolic  equation,  Kolmogorov  spectra 


1 9.  ABSTRACT  (Continue  on  reverse  if  necess^  and  identify  by  block  number) 

Finescale  temperature  features  in  the  oceamc  mixod  layer  (OML)  modify  acoustic  propagation  by  perturbing  the 
index  of  refraction.  A  thermistor  chain  measures  temperature  from  the  surface  to  250  m  with  2  m  vertical  resolution 
and  1.85  m  horizontal  resolution.  Thermistor  data  is  converted  to  high  resolution  sound  velocity  profiles  for  input  to 
an  implicit  flnite  difference  parabolic  equation  (IFD-PE)  model.  The  control  case  is  a  surface  duct  with  a  horizontally 
averaged  sound  velocity  profile.  The  IFD-PE  model  is  run  with  monofrequency  sources  from  500  Hz  to  10  kHz.  The 
acoustic  fluctuations  and  average  acoustic  pressure  are  cennputed  each  20  km  to  a  range  of  100  km.  Acoustic  propa¬ 
gation  through  the  temperature  varying  ON^  is  contrasted  with  acoustic  propagation  through  the  temperature  invari¬ 
ant  OML.  Finescale  temperature  randomly  alters  the  average  acoustic  pressure  by  up  to  a  factor  of  two  and  the  acous¬ 
tic  fluctuations  by  up  to  a  factor  of  five. 


I  [•]  C  7AVJ I  Vj I  m  ^  I : '  TH I 


[3  UNCLASSIFIEDAJNLIMITED  □SAMEASRPT.  □DTlCUSERSl  UNCLASSIFIED 


DD  FORM  1473, 84  MAR  83  APR  sdilion  nwy  be  UMd  unfl  extMUSted  SECURITV  Cl 


83  APR  edilion  nwy  be  ueed  unll  extiausted 
All  o8ier  edWone  are  ofaeotete 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 
LT^CLASSIFIED 


Approved  for  public  release;  distribution  is  unlimited 


COMPUTER  SIMULATION  OF  ACOUSTIC  FLUCTUATIONS 
DUE  TO  FINESCALE  TEMPERATURE  PERTURBATIONS 
MEASURED  BY  THERMISTOR  CHAIN 


by 

Gregory  D.  Crabtree 
Lieutenant,  United  States  Navy 
B.S.,  United  States  Naval  Academy,  1984 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCffiNCE  IN  PHYSICAL  OCEANOGRAPHY 


Author: 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
MARCH  1992 

/  r 

_ r-r'  /V _ _ 

wegfkyD.  Crabtree 


Approved  By: 


Curtis  A.  Collins,  Chairman, 
Department  of  Oceanography 


ABSTRACT 


Finescale  temperature  features  in  the  oceanic  mixed  layer  (OML)  modify  acoustic 
propagation  by  perturbing  the  index  of  refraction.  A  thermistor  chain  measures  temperature 
from  the  surface  to  250  m  with  2  m  vertical  resolution  and  1.85  m  horizontal  resolution. 
Thermistor  data  is  converted  to  high  resolution  sound  velocity  profiles  for  input  to  an  im¬ 
plicit  finite  difference  parabolic  equation  (IFD-PE)  model.  The  control  case  is  a  surface 
duct  with  a  horizontally  averaged  sound  velocity  profile.  The  IFD-PE  model  is  run  with 
monofrequency  sources  from  500  Hz  to  10  kHz.  The  acoustic  fluctuations  and  average 
acoustic  pressure  are  computed  each  20  km  to  a  range  of  100  km.  Acoustic  propagation 
through  the  temperature  varying  OML  is  contrast  with  acoustic  propagation  through  the 
temperature  invariant  OML.  Finescale  temperature  randomly  alters  the  average  acoustic 
pressiue  by  up  to  a  factor  of  two  and  the  acoustic  fluctuations  by  up  to  a  factor  of  five. 
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I.  INTRODUCTION 


Most  acoustic  predictions  are  based  on  strong  assumptions  about  the  oceanic  mixed  layer  (OML)  which 
are  accurate  only  to  a  first  degree.  One  assumption  is  OML  depth  does  not  have  horizontal  variation.  In  fact, 
internal  waves  prcqMigate  along  the  OML-thermocline  interface  and  change  its  dq>th  by  typically  S-lOm. 
Another  assumption  is  the  OML  is  both  isohaline  and  isothermaL  leading  to  a  deterministic  sound  speed 
profile  with  variation  in  depth  dependent  only  on  pressure.  In  reality,  as  cool  underlying  water  is  mixed 
upward  and  warm  surface  water  is  mixed  downward,  there  is  temperature  variation  throughout  the  OML.  The 
temperature  variance  causes  corresponding  variance  in  the  sound  speed  profile.  The  objective  of  this  work  is 
to  detomine  how  good  these  assumptions  are  for  medium  hequencies  jntrpagating  long  distances. 

A  thermistor  chain  data  set  collected  in  October  1989  in  the  Norwegian  Current  system  is  discussed  and 
analyzed  in  Cluster  11.  Thermistor  chain  Hata  provides  a  two  dimensional  profile  from  the  surface  down  to 
2S0  m  with  2  m  resolution  in  the  vertical  and  1.8S  m  resolution  in  the  horizontal.  The  data  is  converted  to  a 
series  of  high  resolution  sound  velocity  profiles  with  spacing  at  every  1.8S  m  in  the  horizontal.  These  profiles 
are  the  input  for  a  computer  simulation  of  acoustic  tran^ission.  Chapter  ni  discusses  acoustic  transmission 
through  a  temperature  varying  OML. 

The  Implicit  Finite  Difference  Parabolic  Equation  (EFD-PE)  acoustic  model  deveicqied  by  Lee  and 
McDaniel  (1987)  is  modified  to  allow  horizontal  variation  at  the  same  distances  as  the  thermistor  chain  data 
set  was  read.  The  IFD-PE  acoustic  model  was  run  at  several  frequencies  from  500  Hz  to  10  kHz.  Two 
simulations  were  run  at  each  frequency:  one  with  a  realistic  OML  and  the  other  with  an  idealized  OML.  In 
the  realistic  case  the  sound  velocity  profile  is  read  into  the  acoustic  model  every  1.85  m.  In  the  idealized  case 
the  temperature  is  averaged  horizontally,  yielding  a  single  sound  velocity  proHle  for  the  entire  event  The 
transmission  loss  with  range  fOT  the  two  cases  is  computed  at  each  frequency.  Average  acoustic  pressure  and 
the  spatial  variance  of  acoustic  pressure  are  shown  in  Chapter  FV.  The  difference  between  sound  prqjagation 
through  a  real  OML  and  through  an  idealized  OML  is  discussed  in  Chapter  V. 

A.  SCALES  OF  MEDIA  SOUND  SPEED  PERTURBATIONS 

The  speed  of  sound  in  the  ocean  is  a  functkxt  of  three  parameters;  temperature,  salinity,  and  pressure. 
An  increase  in  any  one  of  these  three  quantities  causes  a  corresponding  increase  in  the  speed  of  sound.  The 
speed  of  sound  in  sea  water  used  in  this  psper  is  that  determined  by  Mackenzie  (1981) 
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c(z)  =  1448.%  +  4.591  T  -  5304  x  10*  T»  +  2.374  x  10*  T» 

+1.340  (S-35)  +  1.630  x  10*  2  +1.675  x  10*  2*  (1.1) 

-1.025  X  10*  T  (S-35)  -  7.139  x  10"  T  z» 

where  T  is  temperature  in  °C.  S  is  salinity  in  practical  salinity  units,  and  z  is  dq>th  in  meters.  Pressure  effect 
depends  almost  exclusively  on  depth  in  the  water  column  and  pressure  effects  due  to  dynamic  motions  are 
not  considered. 

The  integral  scale  of  convective  mixing  in  the  mixed  layer  is  ^proximately  the  dq}th  of  the  OML. 
Overturning  water  is  an  unstable  {hxx^s  which  creates  turbulence.  There  is  a  cascade  of  energy  from  the 
integral  scale  down  to  1  cycle  per  meto'  (cpm)  as  large  vortices  dissipate  into  increasingly  smaller  vortices. 
Dissipation  dominates  at  distances  less  than  this,  causing  the  cascade  of  temperature  density  spectra  to  decay 
less  rapidly  (Flattd.  1979.  p.  15).  The  region  between  1  m  and  100  m  is  Onestructure  scale  and  is  the  regime 
within  the  purview  of  this  paper. 

Individual  finestructure  features  in  the  ocean  are  large  enough  to  be  accurately  measured  with  modem 
instruments,  but  too  small  and  too  brief  to  be  included  in  ocean  basin  models  and  nowcasts.  Rather  than  model 
isolated  features  indq)endently  it  is  more  practical  to  describe  them  by  wavenumber  spectra.  Temperature 
perturbations  collected  in  this  thermistor  chain  data  set  are  for  this  reason  converted  to  temperature  variance 
spectra.  The  wavenumber  spectra  of  temperature  variance  for  this  analysis  range  between  0.27  m  '  (1  cycle 
every  3.7  m)  and  0.004  m"'  (1  cycle  every  237  m). 

B.  PREVIOUS  WORK 

Underwato'  sound  propagation  has  been  studied  in  earnest  ever  since  World  War  II.  The  earliest 
transmission  loss  models  were  horizontally  homogeneous.  They  relied  on  a  single  sound  velocity  profile 
(SVP)  assumed  to  be  the  same  at  all  locations.  These  models  were  unable  to  accurately  predict  transmission 
loss  across  fronts  and  eddies.  Models  based  on  ray  theory  and  (later)  the  parabolic  equation  ^proximation 
enabled  acousticians  to  model  range  dq}endent  transmission  loss.  Range  indqrendent  and  range  dependent 
models  have  been  relied  upon  to  provide  the  general  features  of  transmission  loss  predictions  right  up  to  the 
present  with  excellent  accuracy.  Their  drawback  is  they  model  a  deterministic  transmission  path  through  a 
very  non-deterministic  medium.  Because  of  the  effects  of  constructive  and  destructive  interference  from 
multipath  arrivals,  slight  variations  in  sound  ^peed  have  a  large  influence  on  the  actual  ranges  of  the  features 
of  transmission  loss  curves.  For  this  reason  acoustic  propagation  is  best  characterized  statistically  for  most 
applications. 
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An  emphasis  on  the  statistics  of  acoustic  prc^jagation  has  generated  considoable  theoretical  and 
experimental  activity.  Theoretical  treatments  have  their  roots  in  scattering  theory  first  pioneered  by  Rytov 
(1937).  He  addressed  scattering  by  smaU  pertuibaticms  of  a  media  valid  for  electromagnetic  pn^gation 
through  the  atmosphere  as  well  as  for  acoustic  propagadon.  Chernov  (1960)  considered  the  case  of  isotropic 
pmurbadons  by  using  the  method  of  small  perturbations  and  ai^lying  an  autoaxrelation  distance  to  the 
inhomogeneities.  His  derivation  relied  on  a  single  tmperature  correlation  coefficier'  rather  than  on  a 
temperature  spectra.  Tatarskii  (1961)  included  a  spectrum  of  scales  of  processes.  Neither  he  nor  Chonov 
specialized  to  ocean  processes  and  they  assumed  media  perturbations  to  be  isotropic.  In  reality  almost  all 
ocean  features  greater  than  1  m  in  length  are  anisotropic,  varying  much  more  in  the  vertical  than  in  the 
horizontal  (Flattf,  1979,  p.  IS). 

Most  of  the  recent  work  on  acoustic  scattering  by  random  media  has  specialized  to  scattering  by  internal 
waves  in  the  deep  sound  channel.  Highlights  are  Munk  and  Zachariason  (1976)  who  used  the  Rytov  method 
(also  called  supereikonal  method)  with  a  gaussian  distribution  of  sound  speed  perturbations  to  address  this 
problem.  They  assumed  the  internal  wave  spectrum  presented  by  Garrett  and  Munk  (1975)  for  the  ocean 
medium.  Flattd,  et  al  (1979)  also  started  with  the  Garrett  •  Munk  internal  wave  model  and  presented 
alternative  acoustic  models  to  describe  acoustic  fluctuations  statistically.  Using  ray  theory,  Ko  (1981)  passed 
multiple  rays  through  a  time  stepped  intenuil  wave  model  within  a  deep  sound  channel.  He  used  the  Monte 
Carlo  method  to  determine  the  statistics  of  acoustic  wave  propagation. 

Although  most  of  the  effort  in  this  field  has  been  expended  on  propagation  through  the  deep  sound 
channel,  there  has  been  an  increased  interest  on  scattoing  of  acoustic  energy  by  Hnestructure  scale 
perturbations.  Wilson  and  Ttqrpeit  (1979)  model  surface  duct  losses  due  to  scattering  by  random  fluctuations 
and  below  duct  ensonification  by  radiation  transport  (i.e.  energy  flux).  Duda,  et  al  (1988)  used  computed 
spectra  of  sound  speed  fluctuations  due  to  microstructure  and  fine  structure  to  determine  the  variance  and  the 
wave  number  spectra  of  acoustic  intensity.  They  investigated  only  the  region  with  small  fiesnel  radius  and 
valid  for  the  Rytov  r^rproximation  which  indicates  that  the  frequency  is  less  than  90  kHz  and  the  range  is  less 
than  1.1  km.  Furthermore,  they  considered  only  a  single  direct  path  propagation  mode. 

Like  many  contributions  to  this  field,  this  thesis  research  is  based  on  a  computer  simulation  of  acoustic 
propagation  to  determine  the  statistics  of  the  propagation.  It  is  unique,  however,  because  it  does  not  rely  on 
theoretical  model  qrectra  of  ocean  features.  It  uses  real  temperature  data.  The  combination  of  an  extremely 
accurate  acoustic  model  and  high  resolution  real  data  provides  a  way  to  directly  compare  acoustic  propagation 
in  a  real  world  horizontally  varying  OML  to  that  of  an  idealized  horizontally  homogeneous  OML,  albeit  with 
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some  simplifying  assumptions.  These  assnnqNions  are  described  in  clugNer  TV.  The  qualitative  results  are 
valid  only  for  the  particular  temperature  perturbadon  distributioo  «k1  simplifying  assumptions  of  this 
simulation.  However,  tiiey  provide  insight  into  the  magnitude  of  the  effects  of  weak  scattering  by 
fuiestructure  qpectia  within  any  surface  duct 
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n.  THERMISTOR  CHAIN  DATA 


A.  descrifuon  of  neat  gin  data  collection 

The  thennistor  chain  data  set  that  this  computer  simulation  is  based  upon  comes  from  the  NEAT  GIN 
ocean  measurement  trial.  NEAT  GIN  data  was  cdlecied  by  Dr.  John  Scott  of  the  United  Kingdom  Admiralty 
Research  Establishment  in  Fall  of  1989.  The  data  set  selected  was  obtained  in  the  Norwegian  Current  system 
on  1 1  October.  It  was  collected  from  1200  to  1800  GMT.  The  sh^  steamed  at  a  nominally  constant  4  knots 
and  sampled  thermistors  every  1.9  seconds  (1.8S2  m)  for  a  total  of  24,000  temperature  soundings.  The  data 
set  traverses  a  44.5  km  straight  track.  Thermistors  were  mouitted  every  2  m  from  the  surface  down  to  250  m. 
This  data  set  provides  an  extremely  high  resolution  two  dimensional  temperature  field. 

B.  DESCRIPTION  OF  THE  MIXED  LAYER  IN  GENERAL 

Turbulent  mixing  forces  the  OML  to  be  almost  isothermal.  Away  from  polar  regions,  temperature 
below  the  OML  typically  decreases  with  depth.  Sound  speed  has  a  maximum  at  the  base  of  the  OML. 
Acoustic  energy  refracts  away  from  the  dq)th  of  the  maximum  sound  speed.  It  is  bent  either  upward  into  the 
duct  or  downward  into  the  deep  ocean.  The  ray  bundles  that  refract  q>waTd  are  reflected  from  the  ocean 
surface  and  continue  to  propagate  within  the  OML.  Since  the  OML  traps  sound,  it  is  also  called  a  surface  duct. 
The  depth  of  the  surface  duct  is  almost  always  the  same  as  the  depth  of  the  OML. 

The  OML  depth  and  temperature  vary  both  temporally  and  spatially.  Figure  1  shows  oscillatory 
feauues  in  the  thermocline  right  below  the  OML.  These  are  caused  by  internal  waves.  Another  characteristic 
of  the  OML  apparent  from  this  figure  is  the  size  of  the  largest  temperature  features.  They  vary  vertically  from 
10  m  to  the  dq>th  of  the  OML  and  vary  horizontally  from  1  km  to  5  km.  Figure  2  makes  the  anisotropy  of 
temperature  features  in  the  OML  even  more  tyjparent  In  this  figure  a  small  segment  of  the  thermistor  chain 
data  is  enlarged  to  show  details.  Two  kilometers  the  thei  .nistor  chain  data  is  displayed  in  Fig  2a).  The  first 
180  m  of  the  same  2  km  segment  is  shown  in  Fig  2b).  They  exhibit  the  anisotropy  of  larger  features  and  the 
near  isotropy  of  the  near  microstructure  scale  features.  Contours  greater  than  40  m  in  depth  traverse  over  1 
km  horizontally.  Those  contours  smaller  than  5  m  in  depth  traverse  less  than  10  m  horizontally. 

Aside  from  horizontal  currents,  most  of  the  kinetic  energy  in  the  world  oceans  is  contained  in  waves. 
This  is  true  everywhere  except  very  near  the  surface  in  the  OML.  Here  turbulence  dominates,  and  motions 
are  regular  in  neither  time  nor  space.  Figures  1  and  2  depict  the  randomness  of  turbulent  kinetic  energy  (TKE) 
in  the  OML. 
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C.  THERMISTOR  CHAIN  DATA  SET 


Shown  below  is  a  contour  {riot  of  temperature  from  the  11  October  1989  thermistor  chain  data 
cdlection.  It  shows  only  the  top  80  m  of  the  entire  M5  km  track.  The  tightly  q)aced  contours  at  the  bottom 


Fig  1.  The  upper  ocean  isotherms  plotted  from  7.9  to  8.2  °C  with  .05  increments. 


of  the  plot  indicate  internal  wave  motions.  They  vary  the  OML  depth  between  70  m  and  83  m.  Temperature 
throughout  the  OML  ranges  0.2  deg  from  8.0  °C  to  8.2  ®C. 

Selection  of  the  11  October  1989  data  set  was  not  based  on  significant  features  or  characteristics 
contained  in  the  data.  Rather,  it  is  suitable  because  it  is  typical  of  mid  and  high  latitude  conditions  through 
much  of  the  year.  The  OML  dqith  remains  relatively  constant  and  temperature  pertuibations  are  mild. 
Modification  of  acoustic  propagation  by  temperature  pertuxbadons  in  this  soun^  channel  is  indicative  of 
I»opagation  through  any  surface  duct. 
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1.  Statistical  Analysis  of  Thermistor  Chain  Data 


a.  HofizontaUy  Averaged  Temperature  Profik 

Most  acoustic  transmission  predictHMis  assume  an  idealized  sound  velocity  profile  based  on 
a  single  temperature  sounding.  They  assume  horizontal  homogeneity  based  on  a  bathythermograph  or 
climatological  data.  This  thermistOT  chain  data  enaUes  cmnparison  of  acoustic  transmission  through  a  real 
OML  and  through  an  idealized  OML.  The  ccMitrol  case  is  an  idealized  OML  generated  by  averaging 
thermistor  chain  values. 

The  horizontally  averaged  temperature  profile  for  the  thermistor  chain  data  collected  on 
October  11. 1989  is  shown  in  Fig  3.  Temperature  perturbations  are  averaged  out,  leaving  a  near  isothermal 


temperature  (°C) 

Fig  3.  Horizontally  averaged  thermistor  chain  data  set.  AD  24,000  profDes  are  averaged  together  at 
each  of  the  74  depths  sampled.  Only  the  top  80  m  is  shown. 
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OML  of  8.1  °C  The  kink  defining  the  depth  of  the  idealized  OML  is  at  70  m.  Below  the  OML  is  a 
thermocline  diat  continues  past  the  250  m  sampled  dq>th.  The  gradient  at  the  top  (rf'  the  theimocline  is  -.05 
°C/m. 

b.  Standard  Deviation  of  Temperature  Versus  Depth 

The  standard  deviation  of  temperature  from  a  horizontal  mean  is  shown  in  Fig  4.  This  figure 


shows  the  standard  deviation  of  temperature  in  the  OML  to  be  about  0.05  °C.  It  jumps  dramatically  below 
the  OML  because  of  internal  waves  propagating  along  the  OML-thennocline  interface.  The  discontinuity  at 
70  m  highlights  the  OML  dqxh.  Temperature  perturbations  above  this  dq)th  ate  characterized  by  small 
deviations  generated  by  turbulent  motions.  The  discontinuity  at  the  depth  of  the  OML  is  slightly  misleading. 
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The  large  temperature  variance  occurs  because  the  thermocline  is  lifted  and  lowered  across  this  depth. 
Horizontal  temperature  changes  in  this  region  are  quite  gradual  because  the  thermocline  is  stratified. 

e.  Temperature  Density  Spectra  at  Selected  Depths 

Temperature  in  the  OML  acts  as  a  passive  contaminant  which  closely  fdlows  turbulent 
velocity  behavior  (Medwin,  1973.  p.  9).  The  Kolmogrmjv  hypothesis  states  that  fully  developed  tuibuloit 
velocity  q)ectra  depcni  only  on  kinematic  shear  viscosity  V  and  the  rate  of  dissipation  of  energy  e.  The 
velocity  spectra  and  by  inference  the  temperature  spectra  O  (k)  have  been  shown  (Batchelor,  1956)  to  be 
related  to  the  wavenumber  of  poturbadons  by 

<I>(A:)  =  Cjk'^ 

where  =  Cj  (V,  £)  is  a  constant  at  each  depth.  The  Kolmogorov  reladonship  is  valid  for  the  inertial 
subregion.  The  temperature  spectra  computed  at  various  d^ths  for  the  thermistor  chain  data  are  shown  in  Fig 

S.  Samples  are  each  237  m  long,  making  the  lowest  resolvable  wavenumber  0.004  m  ^ .  The  Nyquist 

wavenumber  is  0.27  m~^  The  temperature  spectra  fall  within  the  inertial  subregion  fw  which  the 

Kolmogorov  reladonship  is  valid.  A  solid  Une  shows  the  theoredcal  power  law.  At  all  three  depdis 

temperature  variance  is  weighted  primarily  to  the  large  scales  and  decays  by  an  invnse  power  law  to  smaller 
scales.  There  are  no  statisdcally  significant  peaks  within  this  data  that  would  indicate  regularly  varying 
phenomena.  As  expected,  most  of  the  temperature  variance  in  internal  waves  is  associated  with  the  largest 
scales.  Temperature  spectra  in  the  OML  mimic  temperature  spectra  in  the  thermocline,  but  with  much  less 
variance. 
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(m°cO 


Fig  5.  Temperature  spectra  within  the  OML  (30  m),  deep  in  the  theimocline  (250  m),  and  shaDow  in 
the  thennodine  (80  m).  A  solid  line  shows  the  theoretical  -5/3  power  law  for  turbulence  spectra. 


in.  ACOUSTIC  PROPAGATION  IN  INHOMOGENEOUS  WATER 


Acoustic  intensity  \p\  decays  in  a  surface  duct  by  lA  because  of  cylindrical  qrreading.  The  low 
frequency  cutoff  is  qrproximated  by  X.^^y=4.7  x  10”^  where  is  the  wavelength  of  the  cutoff 
frequency  in  feet  and  A  is  the  depth  the  OML  in  feet  (Urick,  1983.  p.  151).  For  this  duct  isS.Omand 
the  cutoff  frequency  is  29S  Hz. 

A.  WAVE  EQUATION  IN  INHOMOGENEOUS  OCEAN 

Before  we  can  detennine  the  effects  of  temperature  perturbations,  it  is  necessary  to  show  the  validity  of 
the  homogeneous  wave  equation  in  a  temperature  varying  surface  duct  First  we  compare  the  influence  of 
pressure  perturbations  on  sound  speed  with  the  influence  of  temperature  perturbations  on  the  speed  of  sound. 
Next  comes  a  derivation  of  the  wave  equation  in  a  temperature  varying  ocean.  The  governing  wave  equation 
in  this  situation  is  not  homogeneous.  The  term  that  makes  the  wave  equation  nonhomogeneous  is  caused  by 
density  differences  in  the  media.  Then  we  show  why  the  variance  of  the  scattered  pressure  due  to  temperature 
change  is  of  greater  importance  than  neglect  of  the  inhomogeneous  term. 

In  the  OML.  the  influence  of  pressure  perturbations  on  the  index  of  refraction  is  much  smaller  than  the 

influence  of  temperature  perturbations  on  the  index  of  refraction.  When  AT  s  0.1  °C. 

J_ 

Co 

perturbations  through  surface  gravity  waves.  When  the  composite  sea  surface  peak  to  trough  Az=3  m. 

]_ 

Co 

assumed  a  function  of  depth  only. 

We  assume  density  of  the  media  to  vary  in  space  and  to  be  frozen  in  time.  This  is  called  “quasi-static” 
because  density  variations  induce  motion,  but  at  a  slower  time  scale  than  the  duration  of  acoustic 
transmissions.  Using  the  method  (rf*  small  perturbations,  we  supoimpose  a  sound  field  over  the  ambient 
density  and  pressure 

p  =  poUy.z)  +pj(jc,y,z) 

P  =  PoU)  +Psix»y>^) 


(dc)  5 

— ^j-Az  is  3.3  X  10  .  For  this  reason  pressure  perturbations  are  not  considered,  and  pressure  is 


(dc)  _4 

(A77  =  3.1  X  10  .  The  greatest  influence  on  the  refraction  index  is  due  to  pressure 
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We  Stan  with  the  equation  of  motion  n^ecting  Coriolis  and  temporarily  disregarding  viscous  effects 


=  -  ^p-pgk  (3.1) 

which  can  be  expanded  as 

(Po  +  Pj)  [d,v  +  0«Vv]  =  -^pQ-Vps- (po+pj)^^ 

Since  p^  «  Pq,  to  the  first  tmler 

Po^r^  =  -  3^  (Pq)  -  ^Ps  -  Po^k 
Subtracting  out  the  hydrostatic  tqjproximation  leaves 

Po3,v  =  -7ps  (3.2) 

The  continuity  equation  states  that 

3^p  +  V  •  (pv)  =  0  (3.3) 

Expanding  this  in  t»ms  of  pg  and  p^  and  retaining  only  rirst  order  terms  leaves 

3,P5  +  V*  (PgV)  =  0  (3.4) 

Differentiating  and  combining  (3.2)  and  (3.4)  yields 

3«Pj  =  0-5) 

In  this  case  the  speed  of  sound  varies  in  space  but  not  in  time  just  as  does  density.  Since  ivessure  is  related  to 
density  by 

dj)  =  cix,y,z)'^d,p 

which  can  be  rewritten 

(Po  +P5)  +  V •  V (po  +Ps)  =  [3, (Po+  P5)  +  V •  V (Pg  +  p^)  ] 

Neglecting  small  terms  leaves 

=  C^[3,(P5)  +  v»Vpg] 

Differentiating  with  respect  to  time  and  rearranging  gives 

^3,,  (Ps)  =  a„  (Pj)  +  0,p)  •  (vp„)  0.6) 

C 

Equations  (3.5)  and  (3.6)  are  now  combined  to  give 


c 

Using  the  relationshv  io  (3*2)  enables  this  to  be  rewritten 

=  ^^Ps~  •'^Ps 

C  Po 

=  ^^Ps-  [V/n(po)  ]  (3.7) 

c 

The  first  two  terms  are  recognizable  as  the  acoustic  wave  equation  in  the  absence  of  temperature 
variability.  The  additional  term  results  from  the  media  doisity  varying  ^latially.  This  equation  is  solved  by 
acoustic  pressure  p^  =  Pq'^Pi  where  Po  is  a  plane  wave 


and  Pj  is  the  scattered  wave.  Chernov  (1960)  shows  diat  the  scattered  wave  is  solved  by 


where 


A  is  the  amplitude  of  the  incident  wave 
k  is  the  wavenumber  of  the  acoustic  transmission 
r  is  the  distance  fiom  a  scattering  element  to  the  receiver  (x,y^)  coordinate 
^  is  the  distance  fiom  source  to  the  scattermg  element,  and 
V  is  the  scattering  volume. 

In  this  case  V  is  the  entire  media  between  source  and  receiver.  Acoustic  fluctuations  due  to  temperature 


perturbations  ate  found  from  the  variance  of  the  scattered  pressure 


— 12 


.This  is 


Here,  N  (r)  describes  the  correlation  distance  of  inbomogeneities  in  the  media.  Both  terms  in  brackets  are 
caused  by  temperature  variance.  The  first  term  in  the  brackets  is  due  to  sound  vdocity  perturbations  and  the 
secondtermin  the  brackets  is  (hie  to  density  poturbations.  Relative  size  calculations  for  the  11  October  data 
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3  1 

show  term  1  »  3.1  x  10  °C~  .  If  the  scattering  angle  6  is  less  than  90  degrees,  tenn  2  is  less  than 

“8  1 

s  9.8  X  10  °C~  and  can  be  disregarded.  Neglect  of  the  second  term  in  the  brackets  of  equation  (3.8) 

is  equivalent  to  neglect  of  [V/n(pQ)  ]  *  Vp^  in  (3.7)  (Ibid,  p.  SO).  This  reduces  the  wave  equation  in  a 
quasi-static  ocean  with  small  temperature  perturbations  to  the  familiar  homogeneous  wave  equation 

\^„iPs)=^^Ps  (3-9) 

C 


B.  PARABOLIC  WAVE  EQUATION 

Separation  of  variables  isolates  the  time  dependence  of  equation  (3.9).  What  is  left  is  the  Helmholtz 
equation.  Dropping  the  subscript  s,  diis  is  written 

V'^p  +  k^p  =  0 

Expressed  in  cylindrical  coordinates,  the  Helmholtz  equation  is 

^r,P  +  7^//^  +  ^zzP  +  (V)  =  0  (3.10) 

whae  kQ  =  (2nf)  / Cq  and  n  (r,  z)  =  Cq/  (c  (r,  z) )  is  *e  local  index  of  refraction.  This  is  an 

elliptic  wave  equation.  Tappert  (1974)  introduced  the  parabolic  equation  approximation  (PE)  to  underwater 
acoustics.  Again,  using  separation  of  variables,  let  the  sound  pressure  field  have  the  form 

pir,z)  =  (3.11) 

Jr 

Here,  1/  (Jr)  contains  the  decay  due  to  cylindrical  spreading.  V  (r)  contains  the  rapidly  varying 
component  of  the  pnessure  field.  It  is  solved  by  the  Hankel  function,  which  is  approximated  in  the  far  field  by 


V  = 


fY' 


A  (r,  z)  is  the  slowly  varying  component  trf'  the  pressure  field.  It  contains  all  of  the  dq>lh  dqtendence. 
A  (r,  z)  is  the  envetopeofsound  pressure  that  the  IFD-re  model  solves.  Applying  the  solution  for  p  (r,  z) 
to  (3.1 1)  disregarding  all  terms  that  diminish  more  rapidly  than  1  /  (  Jr)  gives 


+  +  3,/  +  (^O'*)  =0 

dividing  by  r  ^  V  yields 

kl(n^-i)A+ + likQd^ + a,/  =  0 

If  it  holds  that  «  2ikQ<d^A  then 

kl(n^-l)A  +  lik^d^A  +  a„i4  -  0  (3.12) 

This  is  the  most  widdy  used  parabdic  wave  equation  in  underwater  acoustics.  This  ^qxoximadon  is 
valid  when  the  backscattered  field  is  negligible  and  ray  bundle  travel  is  predominantly  horizontal.  This 
condition  is  easily  met  in  the  OML.  The  PE  accurately  governs  acoustic  prt^tagation  in  a  surface  duct 
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IV.  ACOUSTIC  SIMULATION 


The  acoustic  propagation  model  developed  by  Lee  and  McDaniel  (1987)  iq)plies  a  finite  difference 
numerical  scheme  to  the  of  Tappert  Hnite  differencing  enables  the  governing  equation  to  be  converted 
from  a  boundary  value  problem  to  an  initial  value  problem  in  which  the  solution  marches  horizontally.  This 
means  that  the  pressure  field  at  range  r  +  Ar  is  determined  entirely  by  the  pressure  field  at  range  r.  Because 
the  model  soluticm  marches,  it  is  feasible  to  nnodify  the  soutkI  velocity  profile  at  the  very  short  ranges  of  the 
thermistor  data  collection.  Lee  and  McDaniel  apply  an  “implicit”  finite  difference  algorithm.  Implicit  finite 
difference  schemes  have  the  advantage  of  unconditional  stability.  The  model  is  referred  to  in  this  psqrer  as  the 
IFD-PE. 

The  IFD-PE  acoustic  transmission  model  was  run  at  several  frequencies  between  SOO  Hz  and  10  kHz  in 
both  a  horizontally  temperature  varying  and  temperature  invariant  ocean  medium.  Since  the  IFD-PE  model 
algorithm  matches  horizontally  it  is  especially  well  suited  to  accepting  horizontal  changes  in  sound  speed. 
The  accuracy  of  finite  difference  techniques  is  tied  to  range  step  size.  Lee  and  McDaniel  (1987,  pp.  341-349) 
devote  an  entire  chapter  to  range  stq)  size  analysis.  They  show  an  optimum  choice  for  range  stq)  in  the  IFD- 
PE  is  1/2  of  the  wavelength  of  the  acoustic  frequency.  This  is  the  range  stq>  used  for  all  simulations  of 
acoustic  propagation  in  this  thesis.  For  500  Hz  it  is  1.48  m  and  for  10  kHz  it  is  7  cm.  The  range  step  is 
modified  every  1.832  m  to  make  the  marched  acoustic  pressure  field  coincident  with  the  input  of  each  SVP. 
An  example  of  range  stepping  for  a  500  Hz  source  is  sketched  below. 


1.852  m 


1.852  m 


steprO 


1.480  m  .372  in 
rl  r2 


1.480  m 


.372  m 
r3  r4 


Fig  6.  Step  size  of  the  IFD-PE  model  for  a  500  Hz  sound  source. 


The  same  range  step  size  is  used  for  both  the  varying  and  non  varying  cases  so  that  the  only  difference  in  the 
sound  field  is  due  entirely  to  perturbations  in  sound  speed  caused  by  tempoature  variations. 
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A.  ASSUMPTIONS  IN  THE  COMPUTER  SIMULATION 

1.  TemporaOy  Frozen 

There  are  two  separate  assumptions  about  the  time  variance  of  this  simulatkm.  The  first  is  in  die 
thermistor  chain  data  collection.  An  actual  two  dimensional  temperature  field  must  be  measured  at  all  grid 
points  simultaneously.  This  data  set  is  assumed  to  represent  a  typical  instantaneous  sntqishot  of  an  OML  with 
45  km  in  length.  The  precise  deterministic  paths  of  ray  bundles  are  different,  but  the  statistics  of  acoustic 
propagation  are  as  if  all  tempnature  readings  had  been  taken  simultaneously. 

A  separate  assumption  is  that  the  watH'  is  static  for  the  duration  of  the  sound  propagation  event 
To  understand  this  assumption  it  is  useful  to  sqiarate  water  motions  into  two  categories.  Large  scale  motions 
like  currents  and  tides  advect  the  entire  transmission  path  and  therefore  do  not  affect  acoustic  prqiagation. 
Small  scale  motions  have  a  couple  of  effects.  They  advect  patches  of  water  of  varying  temperature,  thus 
continuously  modifying  the  sound  speed  at  each  location.  They  also  impart  relative  speed  to  acoustic 
transmissions.  The  sound  speed  in  the  ocean  is  about  1500  m/s  and  the  speed  of  small  scale  motions  is 
typically  less  than  10  cm/s.  Although  they  are  five  mders  of  magnitude  smaller  than  the  speed  of  sound,  they 
impart  enough  velocity  to  the  acoustic  energy  to  create  a  doppler  shift  Material  velocity  variations  usually 
have  a  smaller  effect  on  acoustic  transmission  than  do  temperature  variations,  but  they  are  not  altogether 
negligible  (Neubert  1970;  Flatt6.  1979,  p.  86).  The  assumption  of  static  water  motion  is  called  “Taylor's 
hypothesis”  (Flattfi,  1979,  p.  101).  It  is  consistent  with  work  done  by  Chernov  (1960,  p.  36)  among  others. 

2.  Artificial  Absorbing  Bottom 

In  ocean  acoustics  there  are  several  deterministic  paths  ray  bundles  may  follow  between  source 
and  receiver.  In  a  tyincal  ocean  environment  sound  energy  that  re&acts  out  of  the  surface  duct  at  sufficiently 
small  angles  relative  to  the  horizontal  plane  is  refracted  downward  by  the  underlying  thermocline.  Below  the 
thermocline  the  vertical  temperature  gradient  becomes  small.  Furthermore,  the  sound  qreed  increases  with 
pressure,  and  dierefore,  the  sound  speed  increases  with  depth.  The  combination  of  &eep  near-isothomal  water 
and  increasing  pressure  refiracts  rays  toward  the  horizontal  If  there  is  sufficient  depth,  sound  energy  is 
refracted  past  horizontal  and  continues  up  to  the  surface.  This  mode  of  propagation  is  called  convergence 
zone  (C2^.  Sound  energy  that  leaves  the  surface  duct  at  steeper  grazing  angles  interact  with  the  bottom.  There 
they  are  partially  reflected,  partially  absorbed,  and  partially  refracted  through  the  ocean  floor.  Energy 
reflected  frcan  and  refracted  through  the  ocean  bottom  returns  to  the  surface.  This  is  called  the  bottom  bounce 
(BB)  mode  of  prc^tagation.  Both  BB  and  CZ  energy  exit  and  then  reenter  the  surface  duct  The  remaining 
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propagation  paths  are  surface  reflected  at  shallow  angles  and  direct  path  (DP).  These  modes  remain 
continuously  in  the  surface  duct 

CZ  and  BB  pit^gation  is  excluded  fnxn  the  solution  by  introducing  an  artificial  absorbing 
bottom.  This  is  accomplished  in  the  IFD-FC  model  with  an  algorithm  that  extends  the  bottom  by  4/3  of  the 
dqrth  of  the  sound  field  with  a  large  attenuation  coefficient  Acoustic  energy  refracted  or  diffracted  from  the 
surface  duct  continues  downward  because  of  the  thermocline.  It  interacts  with  and  is  attenuated  by  this 
bottom.  This  algorithm  eliminates  all  but  surface  duct  propagation  from  the  solution. 

Inuoduction  of  an  artificially  absorbing  bottom  has  one  drawback  and  three  ovoriding  positive 
aspects.  The  drawback  is  the  artificiality  introduced.  This  is  only  a  minor  difficulty.  Surface  duct  propagaticMi 
is  easily  resolvable  from  BB/CZ  in  the  real  world.  The  primary  advantage  is  that  BB  and  CZ  propagation  do 
not  clutter  the  results.  Since  they  contribute  to  the  ensrmification  of  the  surface  duct  at  intermittent  intervals, 
they  would  be  a  source  of  noise  for  this  simulation.  Another  advantage  is  the  reduction  in  computer  time 
attained  by  modeling  only  the  upper  250  m.  A  16  fold  increase  in  depth  would  result  in  a  16  fold  increase  in 
computer  time.  The  other  advantage  to  a  completely  absorbing  bottom  is  the  difficulty  it  eliminates  by 
determining  how  to  most  realistically  model  the  bottom.  The  coefficient  of  reflection  is  highly  dependent  on 
bottom  type  and  density. 

3.  Gaussian  Starting  Field 

The  PE  model  is  valid  only  in  the  far  field  (r  »  —  )  .  It  needs  a  starting  field  well  removed  from 

ZTC 

the  source.  A  starting  field  can  be  generated  by  a  ray  tracing  model  or  a  normal  mode  solution  to  the 
nonhomogeneous  wave  equation.  A  much  simi^er  initial  field  suitable  fix’  our  purposes,  however,  is  derived 
by  Tappert  (197/).  The  starting  field  consists  of  a  Gaussian  distribution  of  pressure  in  the  vertical.  The 
Gaussian  distribution  is  constrained  by  an  asymptotic  approximation  to  the  solution  of  the  nonhomogeneous 
wave  equation  for  a  point  source.  The  error  induced  by  this  approximation  to  a  point  source  is  negligible  at 
distances  far  from  the  source. 

4.  Isohaline 

As  i»eviously  stated  the  qieed  of  sound  in  water  is  controlled  by  temperature,  salinity,  and 
pressure.  Thermistor  chain  data  provides  the  tempertture  field.  Chemov  (1960)  shows  that  the  effect  of 
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salinity  peruiibadons  cm  sound  speed  is  typically  one  nder  of  magnitude  less  than  that  of  temperature 
pertuitations.  Therefne  we  neglect  the  influence  of  the  salinity  field. 

5.  Other  Scatterers 

There  are  several  other  scatterers  of  acoustic  energy  within  the  OML  not  addressed  in  this 
simulation.  Surface  gravity  waves  and  chop  reflect  sound  energy  out  of  the  OML  and  smear  the  frequency 
transmitted.  Bubbles  attenuate  sound  and  modify  the  local  speed  of  sound.  Material  movement  of  the  medium 
creates  a  doi^ler  shift  that  also  blurs  the  frequency.  Even  biologies  scatter  sound.  These  scatterers  probably 
have  a  greater  influmce  on  sound  qreed  fluctuations  than  do  temperature  perturbations  for  medium 
frequencies.  Medwin  (1973)  found  that  the  predominant  cause  of  acoustic  i^iase  fluctuations  at  24.4  kHz  and 
9S.6  kHz  is  bubble  activity  rather  than  temperature  perturbations.  In  measurements  at  sea,  Nichols  and  Young 
(1968)  found  that  at  270  Hz  there  were  peak  fluctuations  at  the  short  range  (2  nm)  and  the  long  range  (700 
nm)  probably  corresponding  to  surface  waves.  They  also  found  that  for  long  range  only  there  were  peak 
fluctuations  of  0.1  hr'  to  0.01  hr'  probably  corresponding  to  internal  waves.  Urick  and  Tulko  (1%9) 
transmitted  fiequencies  between  7S0  Hz  and  ISOO  Hz  over  24  nm  in  the  deep  sea  to  determine  the  spatial 
correlation  for  bottom  mounted  receiver  separations  up  to  300  ft,  and  they  found  Jiat  the  correlation  was 
unaffected  by  wind  speed  and  sea  roughness,  even  though  the  transmission  path  was  surface  reflected.  Urick 
(1973)  concluded  that  the  major  factors  causing  decorrelation  of  vertical  coherence  from  a  near  surface  source 
are  i)  interference  of  the  direct  path  with  the  reflected  path,  ii)  temperature  and  salinity  perturbations,  and  iii) 
the  effect  of  multiple  surface  and  bottom  reflections. 

B.  STATISTICS  OF  ACOUSTIC  PROPAGATION  THROUGH  OML 

There  are  many  diffoent  ways  to  characterize  the  sound  field.  Units  of  decibels  (db),  acoustic  intensity, 
or  absolute  acoustic  pressure  are  all  equivalent  Because  transmission  loss  generally  vacillates  several  orders 
of  magnitude,  acousticians  often  translate  absolute  acoustic  pressure  P.  to  a  logarithmic  scale  and  use  db.  The 
transmission  loss  in  db  experienced  at  any  (range,  depth)  coevdinate  is  thereby  computed  as  TL=201og(P.). 
In  this  discussion  acoustic  energy  is  everywhere  referenced  to  a  source  of  1  sound  pressure  unit  Absolute 
acoustic  inessure,  therefore,  is  dimensionless.Transmission  loss  curves  are  generally  plotted  at  a  single 
receiver  depth  as  a  function  of  range.  They  display  the  slowly  varying  componem  of  the  sdution  to  the  wave 

I  ikf  ~  cftf) 

equation.  The  rapidly  varying  component  e  is  separated  Grom  the  transmission  loss  curve.  Figure 
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7  is  a  curve  of  this  type.  It  shows  the  decay  of  sound  energy  due  to  cylindrical  spreading  in  the  surface  duct. 
The  sound  pressure  fluctuates  about  a  mean  value  that  is  a  function  of  range.  In  this  case  the  source  and 
receiver  depths  are  both  35  m.  The  curve  has  limited  value  because  it  fails  to  describe  what  the  acoustic 
pressure  is  at  any  other  receiver  depth.  Also,  it  describes  acoustic  fluctuations  inadequately. 

A  more  elaborate  representation  of  the  acoustic  pressure  field  is  given  by  statistical  analysis.  We  need 
both  average  acoustic  pressure  and  the  spatial  variation  of  acoustic  pressure.  Acoustic  pressure  was  recorded 
at  each  2  m  by  2  m  gridpoint  for  400  m  at  each  sample  range  (20  km,  40  km,  60  km,  80  km,  and  100  km). 
Acoustic  pressure  samples  are  only  taken  at  isolated  ranges  rather  than  averaged  together  over  the  whole 
range  of  the  transmission  because  acoustic  pressure  depends  on  range.  That  is,  acoustic  propagation  is  not  a 
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Fig  7.  Normalized  sound  pressure  as  a  function  of  range  for  a  500  Hz  source.  In  both  cases  the 
acoustic  pressure  is  determined  for  a  35  m  source  and  receivers  at  35  m.  The  first  order  decay  of 
acoustic  intensity  with  range  is  due  to  cylindrical  spreading. 


stationaiy  process.  For  purposes  of  comparison,  howev^,  it  is  assumed  to  be  stationary  for  a  400  m  stretdi 
of  ocean. 

The  horizontally  averaged  acoustic  pressure  as  a  function  of  dqxh  is  shown  in  Fig  8  through  Hg  13. 
Each  figure  coneqxMids  to  a  mono&equency  source  from  SOO  Hz  to  10  kHz.  The  source  was  located  at  35  m 
dq>th.  The  ]dots  di^lay  average  absolute  acoustic  pressure  at  range  20  km  and  100  km.  Average  acoustic 
pressure  in  the  idealized  surface  dua  is  shown  with  a  solid  line,  and  in  the  temperature  varying  surface  duct 
is  shown  as  a  dashed  line.  The  figures  show  the  strong  depth  dqiendence  of  acoustic  pressure.  The  maxima 
and  minima  demonstrate  the  depth  dependent  normal  modes.  Higher  frequencies  have  more  modes.  At  some 
frequency-depth-range  combinations,  sound  pressure  is  higher  in  a  temperature  invariant  OML.  At  other 
frequency-depth-range  combinations,  the  situation  is  reversed. 
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Fig  8.  Depth  dependence  of  average  aconstic  pressure  for  a  500  Hz  source  at 
range  a)  20  km  and  b)  100  km.  The  source  is  at  35  m  depth  in  the  surface  duct 
Solid  lines  show  average  aconstic  pressure  in  a  temperature  invariant  OML. 
Dashed  lines  show  average  acoustic  pressure  in  a  temperature  varying  OML. 
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Fig  10.  Depth  dependence  of  average  acoustic  pressure  for  a  2  kHz  source  at 
range  a)  20  km  b)  100  km.  The  source  is  at  35  m  depth  in  the  surface  dnct 
Solid  lines  show  average  acoustic  pressure  in  a  temperature  invariant  OML. 
Dashed  lines  show  average  acoustic  pressure  in  a  temperature  varying  OML. 
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Fig  11.  Dq>tii  dependence  of  average  acoastk  pressure  for  a  3J  kHz  source  at 
range  a)  20  km  lud  b)  100  tan.  The  source  is  at  35  m  depth  in  the  surface  dnct 
Solid  lines  show  average  acoustic  pressure  in  a  temperature  invariant  OML. 
Dadied  lines  show  average  acoosoc  pressure  in  a  temperature  varying  OML. 
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Fig  12.  Depth  dependence  of  average  acoustic  pressure  for  a  5  kHz  source  at 
range  a)  20  km  and  b)  100  km.  The  source  is  at  35  m  depth  in  the  surface  duct 
Solid  lines  show  average  acoustic  pressure  in  a  temperature  invariant  OML. 
Dashed  lines  show  average  acoustic  pressure  in  a  temperature  varying  OML. 
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Fig  13.  Depth  dependence  of  avoragc  acoustic  pressure  for  a  10  kHz  sonrce  at 
range  a)  20  km  and  b)  100  km.  The  source  is  at  35  m  depth  in  flie  surface  duct 


snow  average  acoosoc  pressure  m  a  temperature  mvanant  unit, 
lines  show  average  acoustic  pressure  in  a  temperature  varying  OML. 


The  average  acoustic  pressure  characterizes  the  sound  field  incompletely.  The  qxatial  deviation  of 
acoustic  pressure  is  also  (rf*  interest  Hgure  7  is  an  example  of  how  the  slowly  varying  component  of  the 
solution  to  the  wave  equation  does  not  vary  very  slowly.  It  is  therefne  necessary  to  perform  a  fast  fourier 
transform  on  the  acoustic  [Hessure  at  each  sample  depth  (2m  to  70m).  As  seen  from  Hg  8  through  Fig  13, 
acoustic  pr(q)agation  is  modified  by  depth  as  well  as  range.  However,  for  purposes  of  comparison,  the  spectra 
of  acoustic  fluctuations  is  assumed  to  be  stationary  in  depth  throughout  the  OML  at  each  range.  Averaging 
the  spectra  over  depth  inoeases  the  statistical  relud)ility  of  the  resultant  acoustic  fluctuation  density  curves. 
Any  error  induced  by  vertically  averaging  acoustic  fluctuation  spectra  is  the  same  for  each  case.  An 
assumption  of  ergodicity  at  isolated  ranges  and  through  the  depth  of  the  OML  enables  comparison  of  sound 
fluctuations  in  an  OML  as  a  function  of  range,  frequency,  and  temperature  variance.  They  can  be  compared 
with  a  single  parameter-the  acoustic  fluctuation  density. 

Figure  14  throu^  Fig  19  show  the  acoustic  fluctuation  density  fw  the  six  acoustic  frequencies  of  this 
simulation.  Spectra  of  acoustic  fluctuations  are  displayed  at  ranges  of  20  km  and  100  km.  Again,  solid  lines 
show  computed  results  through  an  idealized  OML,  and  dashed  lines  show  computed  results  through  a  realistic 
OML.  The  spectra  describe  how,  for  a  fixed  instant,  the  slowly  varying  component  of  an  acoustic 
transmission  varies  horizontally. 

Comparison  of  the  acoustic  fluctuation  density  shows  how  acoustic  propagation  is  impacted  by 
horizontal  temperature  variation.  At  frequencies  below  2  kHz,  spectra  of  acoustic  fluctuations  in  a 
temperature  varying  OML  approximate  the  spectra  of  acoustic  fluctuations  in  a  temperature  invariant  OML. 
At  frequencies  above  2  kHz,  the  same  peaks  are  apparent  in  the  spectra,  but  with  much  different  amplitudes. 
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Fig  14.  Energy  density  spectra  of  aconstk  flnctnations  due  to  500  Hz  source  at 
range  a)  20  km,  b)  100  km. 


30 


Fig  15.  Energy  density  spectra 
range  a)  20  km,  b)  100  km. 
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Fig  16.  Energy  density  spectra  of  acoustic  (Inctnatioas  dne  to  2  kHz  source  at 
range  a)  20  km,  b)  100  Im. 
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Fig  17.  Energy  density  q)ectra  of  acoustic  fluctuations  due  to  3.5  kHz  source  at 
range  a)  20  km,  b)100  km. 
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V.  DISCUSSION  OF  RESULTS 


The  sound  field  in  this  computer  simulation  is  a  fwiction  of  sevend  parameters.  It  depends  on  acoustic 
frequency,  range,  and  finescale  structure,  and  source  dqtth.  To  compare  the  relative  importance  of  each 
parameter,  it  is  necessary  to  compute  the  overall  average  acoustic  inessure  and  variance  of  acoustic  pressure 
at  each  frequency  and  range.  The  average  sound  pressure  is  simply  the  mean  of  the  depth  dependent  average 
acoustic  pressure.  The  variance  of  acoustic  fluctuations  is  found  by  integrating  the  spectra  of  acoustic 
fluctuations  over  all  wavenumbers.  This  is  analogous  to  integrating  the  transform  of  a  tinte  series  to  obtain  a 
total  variance.  The  variance  of  acoustic  fluctuations  is  th^fore  depth  averaged  and  wavenumber  integrated. 
It  represents  a  total  spatial  deviation  from  the  mean  acoustic  transmission  loss  through  the  OML. 

A.  Range  Dependence  of  Acoustic  Propagation  in  a  Surface  Duct 

Range  dependence  and  frequency  dependence  of  acoustic  propagation  through  a  surface  duct  are  well 
understood.  To  a  first  approximation,  acoustic  transmissions  decay  with  range  due  to  cylindrical  spreading  in 

a  surface  duct  by  1  /  (  J?)  where  r  is  the  range.  The  average  acoustic  pressure  at  5  ranges  are  plotted  in  Fig 

20.  They  are  compared  with  a  solid  line  showing  decay  with  range  of  1  /  ( Jr) .  More  rapid  decay  is  due  to 
auenuation,  scattering,  and  leakage  out  of  the  surface  duct. 


Fig  20.  Average  acoustic  pressure  versus  range  fai  a  temperature  varying  OML  for  a  1  kHz  source  at 
35  m  depth.  *x*  data  points  show  computed  average  acoustic  pressure.  The  solid  line  shows  the 
theoretical  rate  of  decay  of  acoustic  energy  with  range  due  to  cylindrical  spreading  in  a  surface  duct. 
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Fig  21.  a)  Variance  of  acoustic  fluctuations  with  range  for  a  1  kHz  source  in  a  temperature 
varying  OML  at  35  m  depth.  In  b)  the  standard  deviation  of  acoustic  fluctuathms  fe  compared  to 
a  s<^  one  showing  the  rate  of  deoiy  of  average  acoustic  energy  due  to  geometrical  spreading  in  a 
surface  duct 


Acoustic  fluctuations  appear  to  decay  witt  range  similar  to  the  cylindrical  spreading  of  average 
transmission  loss.  Figure  21a)  has  the  acoustic  fluctuation  variance  plotted  for  a  1  kHz  source  at  35  m  depth 
in  a  temperature  varying  OML.  Variance  is  computed  every  20  km  from  20  km  to  100  km.  Because  it  has  the 
same  units,  a  better  parameter  to  compare  with  average  transmission  loss  is  die  standard  deviation  of  acoustic 
fluctuadons.  Hgure  21b)  shows  the  standard  deviation  of  acousdc  fluctuadons.  It  is  compared  to  a  solid  line 

showing  decay  of  1/  ijr)  .The  solid  line  passes  through  the  first  value  of  standard  deviation  at  20  km.  This 
figure  demonstrates  that  acousdc  fluctuadons  dimiiush  less  with  range  than  average  transmission  loss  does. 
Similar  comparisons  show  that  the  ratio  of  standard  deviation  of  sound  pressure  to  average  sound  pressure 
increases  with  range  at  all  frequencies. 

B.  Frequency  Dependence  of  Acoustic  Propagation  through  a  Surface  Duct 

Variance  of  acousdc  fluctuadons  in  this  surface  duct  shows  a  predictable  dqiendence  on  frequency.  It 
has  a  maximum  value  that  is  a  function  of  range  as  well  as  frequency.  Figure  22  shows  frequency  dependence 
at  20  km  and  at  100  km.  The  source  is  at  35  m  depth  and  sound  propagates  through  a  temperature  varying 
OML.  At  20  km,  variance  of  acoustic  fluctuations  is  the  highest  at  3.5  kHz.  It  is  less  at  frequencies  both  higher 
and  lower.  The  greatest  variance  at  100  km  is  at  1  kHz.  This  trend  is  predictable  because  it  is  very  similar  to 
the  dependence  of  average  transmission  loss  on  frequency  in  a  surface  duct.  Acoustic  energy  is  leaked  out  of 
the  duct  at  relatively  low  frequencies  near  the  duct  cutoff  frequency.  Acousdc  energy  at  higher  frequencies  is 
lost  by  attenuation.  Since  acoustic  fluctuations  are  driven  by  average  acoustic  energy,  it  is  anticipated  that 
variance  of  sound  transmission  adh^e  to  similar  frequency  dq)endence. 
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Fig  22.  Variance  of  acoustic  fluctuations  vmus  freqnencj  at  range  a)  20  km  and  b) 

100  km.  The  sound  source  is  at  35  m  and  sound  transmits  through  a  temperature 
varying  OML. 

C.  Dependence  of  Acoustic  Propagation  on  Flnescale  Temperature  Perturbations 

The  stated  objective  of  this  computer  simulation  has  been  to  isolate  all  parameters  affecting  sound 
transmission  to  compare  propagation  through  a  real,  temperature  varying  OML  with  imqiagation  through  a 
temperature  invariant,  idealized  OML.  Figure  8  duough  Fig  19  in  Ch^iter  IV  indicate  that  acoustic 
transmission  is  markedly  different  through  the  two  mixed  layers.  The  qiectra  of  acoustic  fluctuations  are 
impacted  the  most  at  frequencies  above  2  kHz.  The  overall  variance  between  the  two  cases  at  emtain 
frequoKies  differs  by  as  much  as  a  factor  of  S.  Howevn,  temperature  spectra  do  not  modify  either  the  average 
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acoustic  pressuie  or  the  acoustic  fluctuation  qwcira  in  a  predictid>le  paaem.  Figure  23  shows  average 
acoustic  pressure,  and  Figure  24  shows  variance  of  acoustic  fluctuations  plotted  against  frequency  at  each 
range  analyzed.  The  data  points  marked  ‘x’  are  variance  caused  by  propagation  of  a  source  at  35  m  dq)ih 
through  a  temperature  varying  envitonmenL  The  data  points  marked  ‘o’  are  die  same  except  transmission  is 
through  a  temperature  invariant  OML.  There  is  no  clear  trend.  The  probaUe  cause  of  this  is  the  strong 
dqwndence  of  sound  on  OML  depth. 

D.  Dependence  of  Acoustic  Propagation  on  Source  Depth 

To  investigate  dqxh  dqiendence,  the  source  depth  of  35  m  was  modified  at  3.5  kHz  and  at  5  kHz.  The 
IFD-PE  model  simulated  acoustic  transmisskm  through  an  idealized  OML  and  a  realistic  OML  at  3.5  kHz 
with  source  located  at  5  m.  Simulations  were  also  run  fcv  a  5  kHz  source  located  at  5  m,  30  m,  and  45  m.  The 
dq)th  dependence  of  average  acoustic  {vessure  is  shown  in  Hg  25,  and  the  depth  dependence  of  acoustic 
fluctuation  variance  is  shown  in  Fig  26.  Holding  all  parameters  cmistant  excqM  dq)th  can  alter  the  variance 
of  acoustic  fluctuations  by  a  factor  of  5.  Translating  the  sound  source  just  5  m  from  35  m  to  30  m  increases 
acoustic  flucoiations  at  60  km  in  a  temperature  invariant  OML  three  fold.  This  is  evidence  that  the  precise 
combination  of  source  depth,  frequency,  and  OML  depth  have  a  greater  influence  on  acoustic  transmission 
than  does  temperature  poturbations.  The  source  depdi,  frequency,  and  OML  depth  combination  activate 
particular  normal  modes,  which  in  turn  raise  and  lower  the  average  acoustic  transmission  and  acoustic 
fluctuation  density.  This  explains  why  temperature  spectra  do  not  modify  the  sound  field  in  a  predictable 
pattern.  The  realistic  OML  depth  varies  with  range,  changing  which  modes  are  activated.  The  change  in 
geometry  of  the  frequency  wavelength  compared  to  the  source  depth  and  OML  depth  with  range  is  a 
stochastic  process. 
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Fig  23.  Depth  averaged  acoustic  pressure  as  a  function  of  frequency  at  range  a) 
20  km,  b)  40  km.  Source  depth  is  35  m.  *o’  data  points  indicate  results  for  the 
ideali^  OML.  ‘x’  data  potaits  are  for  the  temperature  varying  OML. 
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Fig  23.  (coot)  I>q)tii  averaged  acoustic  preasure  as  a  ftaactiou  offt^queiicy  at 
range  c)  M  km  and  d)  80  km.  Source  deptk  is35  m.  *o’  data  points  indicate 
res^  for  the  idealized  OML.  *x’  data  points  are  for  the  temperature  varying 
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Fig  23.  (cont)  Depth  averaged  acoustic  pressure  as  a  functhm  of  frequency  at 
range  e)  100  km.  Source  depth  is  35  m.  ‘o*  data  pants  indicate  results  for  the 
idealized  OML.  ‘x’  data  points  are  for  the  temperature  varying  OML. 
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Fig  24.  Variance  of  aconstk  flnctnatiou  as  a  function  of  fircqnency  at  range  a) 
20  km,  b)  40  km.  Source  dinth  is  35  m.  'o*  data  points  indicate  results  for  die 
ideali^  OML.  ‘z’  data  pomts  are  for  the  temperature  varying  OML. 
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Fig  24.  (cont)  Variance  of  acoustic  fluctuations  as  a  function  of  frequency  at 
range  c)  60  km  and  d)  80  km.  Source  depth  is  35  m.  *o’  data  points  indicate 
results  for  the  idealized  OML.  ‘x*  data  points  are  for  the  temperature  varying 
OML. 


dqidi(ni) 

Fig  25.  Dependence  of  average  acoustic  pimure  on  source  depth  for  a  5  kHz 
source,  ‘x*  indicates  propagation  through  a  temperature  varyfaig  OML.  *o* 
indicates  propagation  thro^  a  temperature  invariant  OML. 
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Fig  24.  Dependence  of  variance  of  acoustic  floctnatioiisoDdqitfa  for  a 
monoftrequency  source  at  a)  kHz,  and  b)  5  kHz.  *x’  indicates  propagation 

throngh  a  temperature  varying  OML.  V  indicates  propagation  throng  a 
temperature  invariant  OML. 


VL  CONCLUSIONS 


A.  HIGHUGHTS 

The  IFD-PE  model  is  a  valuable  tool  to  contrast  acoustic  transmission  duough  a  temperature  varying 
OML  with  acoustic  pnqngation  through  a  temperature  invariant  OML.  A  computer  experiment  was  designed 
using  the  transmission  model  that  iscdated  all  parameters  capable  of  influencing  sound  propagation.  A 
statistical  analysis  of  the  sound  field  enabled  a  direct  comparison  of  the  impact  of  temperature  perturbations 
on  acoustic  int^gation.  Depth  icpeadeat  average  acoustic  pressure  was  found  to  vary  by  up  to  a  factor  of 
two  when  crmtrasted  for  the  two  cases.  Spatial  variance  of  acoustic  fluctuations  was  found  to  vary  by  iq>  to  a 
factor  of  five.  Acoustic  fluctuations  are  closely  related  to  average  sound  pressure.  The  standard  deviation  of 
spatial  fluctuations  was  found  to  be  on  the  order  of  10  percent  of  the  depth  averaged  acoustic  pressure.  It  is 
lower  for  low  frequencies  and  much  higher  fOT  high  ftequencies.  Additional  simulations  with  changes  in 
source  depth  suggest  that  the  impact  of  real  temperature  profiles  on  the  sound  field  is  indirect  It  randomly 
modifies  the  normal  modes  that  can  propagate.  This  thesis  proves  that  fw  medium  acoustic  firequencies 
propagating  through  a  typical  OML,  an  assumption  of  isothennal  water  and  constant  dqpth  are  valid. 

B.  SUGGESTIONS  FOR  FURTHER  WORK 

Additional  effort  is  needed  to  verify  the  statistics  of  an  ensonified  surface  duct  The  approximate  values 
of  a  factor  of  two  for  average  acoustic  pressure  and  r^o  of  one  to  ten  for  standard  deviation  to  average 
acoustic  pressure  may  be  OML  dqpendent  Also,  no  consideration  was  given  to  the  horizontal  dqrendence  of 
temperature  fluctuations.  Each  simulation  started  with  the  same  sound  velocity  profile.  If  sound  propagation 
is  not  modified  by  horizontal  changes  in  source  location,  then  temperature  spectra  akme  are  responsible  for 
altering  the  sound  field.  If  hcxizonial  changes  in  source  location  alter  the  sound  field,  then  precise  temperature 
features  are  re^nsible. 


48 


LIST  OF  REFERENCES 


Batchelor,  G.  K.,  et  al,  “Small-scale  variation  of  convected  quantities  like  temperature  in  tuitwlent  fluid” 
(Pans  1  and  2)./.  Fluid Mech..  34.  pp.  443-448. 1959. 

Chernov,  L.  A.,  Wave  Propagation  in  a  Random  Medium  (translated  by  R.  A.  Silvennan),:  Dover 
Publications.  New  Yrak.  1960. 

Duda,  T.  F.,  et  al,  “Modelling  meter-scale  acoustic  intenaty  fluctuations  from  oceanic  fine  structure  and 
microstructure,’’  J.  Geophysical  Research,  93.  pp.  5130-5142, 1988. 

Flattf,  S.  M.  ed.,  Sound  Transmission  Through  a  Fluctuating  Ocean,  Cambridge  University.  New  York, 
1979. 

Lee,  D.  and  S.  T.  McDaniel,  Ocean  Acoustic  Propagation  by  Finite  Difference  Methods,  Peigamon  Press, 
Oxford,  1987. 

Mackenzie,  K.  V.,  “Equation  of  sound  speed  in  the  ocean,”  J.  acoust  Soc.  Am.,  70.  p.  808, 1981. 

Medwin,  H.,  Predicting  Sound  Phase  and  Amplitude  Fluctuations  due  to  Microstructure  in  the  Upper 
Ocean,  unpublished  NPS  technical  report  NPS-61Md73111  A,  1973. 

Neubert,  J.  A.  and  J.  L.  Lumley,  “Derivation  of  the  stochastic  helmholtz  equation  fiv  sound  propagation  in  a 
turbulent  fluid”,  J.  acoust  Soc.  Am.,  48,  pp.  1212-1218, 1970. 

Tappea,  F.  D.,  ‘Tarabolic  equation  method  in  undowater  acoustics,”  J.  acoust.  Soc.  Am.,  55,  p.  534(A), 
1974. 

Tappen,  F.  D.,  “The  parabolic  equation  method”,  Wave  propagation  and  Underwater  Acoustics  (Eds.  J.  B. 
Keller  and  J.  S.  Papodakis);  Lecture  Notes  in  Physics,  70,  Springer,  Berlin,  1977. 

Tatarskii,  V.  I.,  The  Effects  of  the  Turbulent  Atmosphere  on  Wave  Propagation  Gsrael  Program  for  Scientific 
Translation),  1971. 

Urick,  R.  J.,  “Measurements  of  the  vertical  coherence  of  the  sound  from  a  near-surface  source  in  the  sea  and 
the  effects  on  an  additive  vertical  array,”  /.  acoust.  Soc.  Am.,  54,  pp.  1 15-120, 1973. 

Urick,  R.  J,  Principles  of  Underwater  Sound,  3rd  Ed,  McGraw-Hill,  New  York,  1983. 

Urick,  R.  J.  and  T.  J.  Tulko,  “Vertical  coherence  of  sound  transmitted  over  a  twenty  four  mile  path,”  J. 
Acoust.  Soc.  Am.,  46,  pp.  1308-1317, 1%9. 


49 


INITIAL  DISTRIBUTIOri  LIST 


Defense  Technical  Infomaation  Center  2 

Cameron  Station 
Alexandria.  VA  22304-6145 

Dudley  Knox  Library  2 

Code  052 

Naval  Postgraduate  School 
Monterey,  CA  93943 

Chairman,  Code  OC  1 

Oceanography  Department 
Naval  Postgr^uate  School 
Monterey,  CA  93943 

Assoc.  Professor  P.C.  Chu  1 

Oceanography  Department 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Assoc.  Professor  A.A.  Atchley  1 

Physics  Department 
Naval  Postgraduate  School 
Monterey,  CA  93943 

Dr.  J.C.  Scott  1 

Admiralty  Research  Establishment 

Southwell,  Portland 

Dorset  DT5  2JS 

United  Kingdom 

Dr.  D.S.  Ko  1 

Institute  for  Naval  Oceanography 
Stennis  Space  Center,  MS  39529-5005 

Oceanographer  of  the  Navy  1 

Naval  Observatory 

34th  and  Massachusetts  Avenue  NW 

Washington,  DC  20390-5000 


50 


Dr.  D.  Lee 

Naval  Underwato'  Systems  Center 
New  London,  CT  06320 

Professor  S.M  Flatt6 

Physics  Department 

University  of  California  -  Santa  Cruz 

Santa  Cruz,  CA  95064 

Professor  H.  Medwin 
Physics  Department 
Naval  Postgraduate  School 
Monterey,  CA  93943 

CDR  J.  Fink 

Commander  Task  Force  Twelve 
Pearl  Harbor,  HI  96601 

LT  G.D.  Crabtree 
Naval  Air  Facility  ASWOC 
Lajes  Field,  Azores  09720 


51 


